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Radio images of supernova remnants in the framework of diffusion model are discussed. The distribution 
profiles of synchrotron radiation intensity for spherical injection source of relativistic electrons are reduced 
at different frequencies. An explanation of the observational data obtained on UTR-2, according to which 
the size of the supernova remnant at decametric waves is larger than the remnant size at high frequencies, 
is given. 

Introduction 

The radio images of the supernova remnants (SNR) obtained with the UTR-2 radio telescope by M.A. 
Sidorchuk and E.A. Abramenkov show the difference between SNR sizes at low frequencies (decametric 
waves) vs. the higher ones [1]. 

For example, the brightness temperature contour map of SNR HB3 at the frequency 20 MHz is shown 
in Fig.l. The temperatures on the isophotes are given in 10 3 K. At the lower corner the UTR-2 directional 
diagram is shown as a dashed circle. The center of SNR is denoted by the white " +". At the frequency 1420 
MHz radio image of SNR corresponds to the size 60 x 80 pc (the region limited by a white circumference) [2J. 
At frequency 20 MHz it corresponds to the size 70 x 95 pc pj. We can see that the SNR size at low frequency 
is larger than the SNR size at high frequency. In this work we propose an explanation of the SNR size 
changing with frequency in the framework of diffusion model. 

Diffusion model 

In the diffusion model we consider the SNR as a spherical source of relativistic electrons. The particles 
are accelerated by the spherical shock wave (S W) front and next further propagate by means of diffusion 
into the environment medium losing their energy due to synchrotron radiation [3]. The basic equation is the 
kinetic equation (KE) for electron distribution function (EDF) N(E,t,f) with a source Q(E,t,r) and with 
diffusion as a mode of the electron propagation 

were A is the Laplace operator, D is the diffusion factor that as we will show bellow can be chosen to have 
no energy dependence. 

The second term in KE describes the synchrotron losses of electron 
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where e is the electron charge, m is the electron mass, c is the light velocity and Wjj. = H 2 /(8ir) denotes the 
energy density of the magnetic field. 

The right part of equation JT]) corresponds to the source of relativistic electrons: 

Q(E,t,r) = Q • N mj {E) ■ S(r) • 6(t), 



36 



YbU ±4 of Contributed rapers 



A. V. Karnausaenko, £,. Yu. aannikova, V. M. Kontorovich 



where S(f) = S(\r\ — R) sets the spherical form of the source, with radius R, Ni n j = E 70 • Q(E max — E) ■ 
Q(E — E m i n ) is the distribution function of injected electrons, 70 = 2, G(...) is the Havyside's step function. 
The solution of KE is the EDF for the spherical source [3], pQ: 

N(E,t,r) = ®° 2 ■ dr'_ • ^ + ^~\ • N inj (E) • R* [ dtt • exp ( ^ ~ f) \ ) , (2) 
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where 
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(r-r 1 ') 2 = (x-x') 2 + (y-y') 2 + (z-2; / ) 2 , x' = i? • sin • cos ip, y' = R ■ sin 9 sin </?, z' = R- cos 9, dfl = 
sin 9d9d<p. 

The physical parameters of the diffusion model are the following: he life time r = l/((3 ■ E) - the time of 
particle existence with energy E, diffusion length l^if = \JD$ ■ t - the distance, which particle passes during 
the life time, and the diffusion velocity is Vdif = Do/l^if. The diffusion length defines the size of the radio 
image of SNR. Knowing the EDF, we can calculate the intensity distribution of the synchrotron radiation / 
[5]. The intensity distribution in the image plane for the sources resolved by the radio telescopes is given by: 

10, t, r) = ^ 3- / dE dz- N(E, t,f)-H ± --- K 5/3 (r,)d V , (3) 
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where H± is the magnetic field projection on the image plane, J dz is an integral along the line of sight and 
A'5/3 is modified Bessel function. The numerical calculations were made with help of the Mathematica 5.1 
package. 



Frequency dependence of synchrotron radiation intensity distribution 

The profiles of synchrotron radiation intensity distribution at different frequencies are given in Fig. 2. 
Though the remnant is not fully symmetrical we choose a sphere with R = 50pc as approach. We neglect 
the motion of the shock front and consider the SNR as an immovable injection source. 

In Fig. 2a r = 3.6 • 10 11 s,/^j = 1.3 • 10 20 cm, Vdif = 3.7 • W 8 cm/s: the diffusion length is higher than 
SNR radius, it almost completely defines the size of the radio image. In Fig. 2b r = 8 • 10 10 sj^if = 
6.3- 10 19 cm, Vdif = 7.8- 10 8 cm/ s, the life time and the diffusion length decrease, and a gap (valley) appears in 
the figure. In Fig. 2c the life time and the diffusion length further decrease r = 5-10 10 s, Idif = 5-10 19 cm, Vdif = 
9 • 10 8 cm/s, and the gap increases. In Fig. 2d with parameters r = 4.2 • 10 10 s,^j/ = 4.6 • 10 19 cm, Vdif = 

1 • 10 9 'cm/s the gap reaches maximum value: the diffusion length becomes less than own size of the system 
and has no influence in its radio image. The approximation of a constant diffusion factor is valid if the 
following equality is carried out: Idif 1 /Idif 2 ~ \Jvi/v\. For the chosen frequencies and corresponding values 
of diffusion lengths {l di fi = 1.3 • 10 20 cm, v x = 20MHz and l di f 2 = 4.6 • 10 19 cm, v 2 = U2QMHz) the equality 
is carried out with a good accuracy 1.3 • 10 20 cm/4.6 • 10 19 cm \J lA20MHz / 20MHz « 2.9 It is necessary 
to note, that velocity of the SNR shock front R(t) = 3 • 10 7 cm/s can be of the same order as the diffusion 
velocity and even more, therefore in general it is necessary to consider the front movement. 



Conclusions 

In the framework of the diffusion model we have explained the observational data obtained by the UTR-2 
according to which the size of the supernova remnants on the low frequencies is is higher than the one at 
high frequencies. It has been shown that in our case the diffusion factor may be considered as a constant. 
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Figure 1. Brightness temperature contour map of SNR HB3 at frequency 20 MHz by M.A. Sidorchuk 
and E.A. Abramenkov [T|. 



Figure 2. Profiles of intensity of synchrotron radiation distribution for spherical source with constant 
diffusion factor for following parameters a) v = 20MHz; b) v = A08MHz; c) v = 1000MHz; d) v = 
1420MHz, I u is expressed in relative units in the same scale. 



Figures are available on YSC home page (http://ysc.kiev.ua/abs/procl4_8.pdf). 
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